A global fit study on the new agegraphic dark energy model 
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We perform a global fit study on the new agegraphic dark energy (NADE) model in a non-flat universe 
by using the MCMC method with the full CMB power spectra data from the WMAP 7-yr observations, the 
SNIa data from Union2.1 sample, BAO data from SDSS DRV and WiggleZ Dark Energy Survey, and the latest 
measurements of Hq from HST. We find that the value of ^^^o is greater than at least at the 3cr confidence levels 
(CLs), which implies that the NADE model distinctly favors an open universe. Besides, our results show that 
the value of the key parameter of NADE model, n = 2.673+[J ;Jj^^+^ ^g, at the \-3o- CLs, where its best-fit 
value is significantly smaller than those obtained in previous works. We find that the reason leading to such a 
change comes from the diff'erent SNIa samples used. Our further test indicates that there is a distinct tension 
between the Union2 sample of SNIa and other observations, and the tension will be relieved once the Union2 
sample is replaced by the Union2. 1 sample. So, the new constraint result of the NADE model obtained in this 
work is more reasonable than before. 

PACS numbers: 95.36.+X, 98.80.Es, 98.80.-k 



I. INTRODUCTION 

Dark energy has become one of the most important research 
areas in cosmology since the cosmic acceleration was discov- 
ered by the observations of type la supernovae (SNIa) [ 1] . The 
cosmological constant, A, introduced by Einstein in 1917, is 
a natural candidate for dark energy, which fits the observa- 
tional data best to date. However, the model of A plus cold 
dark matter, known as the ACDM model, always suflTers from 
the two well-known theoretical difl^iculties, namely, the fine- 
tuning and cosmic coincidence problems [2]. Thus, dynami- 
cal dark energy models become popular, because they may al- 
leviate the theoretical challenges faced by the ACDM model. 
For a recent review of dark energy, see Ref. [3]. 

Among the many dynamical dark energy models, a class 
with feature of quantum gravity, embodying holographic prin- 
ciple, looks very special and attractive. Such class of mod- 
els, usually called "holographic dark energy", has been in- 
vestigated in detail [4] and constrained by the observational 
data [5]. In this paper, we focus on the "new agegraphic dark 
energy" (NADE) model [6] in this class, which takes into ac- 
count the uncertainty relation of quantum mechanics together 
with the gravitational eflTect in general relativity. The reason 
that the NADE model is fairly attractive is also owing to the 
fact that this model has the same number of parameters to the 
ACDM model, less than other dynamical dark energy mod- 
els. The NADE model has been proven to fit the data well. 
However, it should be noticed that the global fit analysis on 
the NADE model has never been done. 

The global fit analysis is very important for the study 
of a cosmological model. To determine the parameters of 
dark energy models, one often has them constrained by us- 
ing some important observational data sets, such as SNIa, 
baryon acoustic oscillation (BAO) and cosmic microwave 
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background (CMB). For the CMB information, to save time 
and power, one often uses the CMB "distance priors" (R, 
I A, z*) whose information is extracted from the temperature 
power spectrum of CMB based on a ACDM scenario. How- 
ever, it should be noted that using the CMB "distance priors" 
to constrain other dark energy models will lead to a circular 
problem, since the values of the CMB "distance priors" de- 
pend on a specific cosmological model, namely, the ACDM 
model. Besides, the "distance priors" do not capture the infor- 
mation concerning the growth of structure probed by the late- 
time ISW eflTect which is sensitive to the properties of dark 
energy. In the past, the NADE model was constrained only by 
the CMB "distance priors" data. In this study, we shall use the 
full CMB data from the 7-yr WMAP observations to explore 
the parameter space of the cosmological model involving the 
NADE. Such an exploration is believed to complete the study 
of the NADE model. 

In the previous studies, the spatial curvature, Qy^o, is usu- 
ally neglected. So, we do not know if a flat universe is still 
favored by current data in the NADE model. Moreover, the 
initial condition of the NADE model [6] is only applicable 
for a flat universe containing only dark energy and pressure- 
less matter. In Ref. [7], a new, updated initial condition is 
proposed, which is applicable for a non-flat universe with var- 
ious components. We shall apply this new initial condition 
in the present work. In such an application, for employing 
the "CosmoMC" package, one may treat the parameter n as 
a derived parameter and use the Newton's iteration algorithm 
to determine its value. In addition, in a global fit analysis, 
the cosmological perturbations of dark energy should also be 
considered. Thanks to the quintessence-like property of the 
NADE, there is no w = -1 -crossing behavior and so the cor- 
responding gravity instability is absent. In the holographic 
dark energy model, the w = -1 -crossing happens, so one has 
to employ the parameterized post-Friedmann (PPF) approach 
to treat the gravity instability of the perturbations in a global 
fit analysis; for details see Ref. [8]. For the NADE model, 
however, we have no need to employ the PPF approach. 

Our paper is organized as follows. In Sec. II, we derive 
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the background and perturbation evolution equations for the 
NADE model in a non-flat universe. In Sec. Ill, we introduce 
the fit method and the observational data, and then give the 
global fit results. Concluding remarks are given in Sec. IV. 
In this work, we assume today's scale factor gq = I, and so 
the redshift z satisfies z = - I; the subscript "0" always 
indicates the present value of the corresponding quantity, and 
the units with c = fi = 1 SLve used. 



II. BRIEF DESCRIPTION OF THE NADE MODEL 

The NADE model is constructed in light of the Karolyhazy 
relation and corresponding energy fluctuations of space-time. 
The energy density of NADE, pde, has the form [6] 



Pde 



^2 



(1) 



where ^ is a numerical parameter, Mpi is the reduced Planck 
mass, and t] is the conformal age of the universe. 



r'dl_ r da 
^ Jo a Jo Ha^' 



(2) 



where a is the scale factor and H = a/a is the Hubble pa- 
rameter. Here a dot denotes the derivative with respect to the 
cosmic time t. 

The background evolution of a non-flat universe is de- 
scribed by the Friedmann equation. 



E 



a = 1, 



(3) 



where is defined as the ratio of the energy density p/ to 
the critical energy density pc = SMp^H^ with / = de, dm, b, 
r and k which denotes the dark energy, dark matter, baryon, 
radiation and spatial curvature, respectively. Note that here 
Pk = -3Mlik/a^. 

The energy conservation equations for the various compo- 
nents take the form 



A + 37^(l-hw0p/ = O, 



(4) 



where w/ is the equation-of- state parameter (EOS) of a spe- 
cific component denoted by subscript /, such as Wdm = wt = 0, 
Wr = 1/3 and Wk = -1/3. For the EOS of NADE, one can ob- 
tain 



, , 2(1 +Z) 

-I + — V^^^^ 

3n 



from Eqs. (1), (2) and (4). Furthermore, using Eqs. (l)-(4), 
one can also easily get the equation of motion of Qj^, a diff'er- 
ential equation. 



For solving the difl'erential equation (6), we adopt the initial 

condition, ^deiZini) = ^^^^^1^ (l -h ^/F(E~i)f at zm = 2000 

with F(z) = Q ^+Qq(i+z) ' giv^i^ ill f^^f- [^1' which helps reduce 
one free parameter for the NADE model. This initial condition 
is an updated version of the original one, flde(Zini) = n^(l -\- 
Zini)~^/^ at Zini = 2000, given in Ref. [9]. The only difl'erence 
between them is that the new one is valid when the radiation 
is considered as a non-ignorable component. Note that the 
initial condition is also applicable in a non-flat universe, since 
the spatial curvature flk is much smaller than fl^ or at z = 
2000. For more details about the derivation and application of 
the new initial condition, see Ref. [7]. 

In our work, for simplicity, we fix D^o = 2.469x 10"^/z"^(l -i- 
0.2271A^eff) according to the WMAP 7-yr observations [10], 
where A^eff = 3.04 is the standard value of the eff'ective number 
of neutrino species and h is the Hubble constant Hq in units of 
lOOkm/s/Mpc. 

It should be noted that the parameters n, ^dmo^ ^bo and 
Q^ko are not independent. One of them can be derived from 
the others. For example, we can treat the parameter n as 
the derived parameter. In this case, once the values of fldmo, 
flbo and flko are given, one can impose an initial value of n 
and then obtain the true value of n by using the Newton it- 
eration algorithm, rik+i = rik - f(nk)/ f(nk), where f(nk) = 
1 - ^dmo - ^bo - ^ro - ^ko " ^de(0)\n, with ^^^(0)^, the ^th 
order numerical solution (at z = 0) of Eq. (6) with the new ini- 
tial condition, and f\nk) is the numerical derivative of f(n) at 
the point n = rik. Of course, one can also treat ^dmo or Q/,o or 
flko as a derived parameter and get its value by using a similar 
method. 

For our global fit analysis, it is convenient to modify the 
"CosmoMC" package [11] if ^ is set as a derived parame- 
ter. Thus, we will use the Newton method mentioned above, 
which is proved to be convergent, to get the true value of n. 
We choose the initial value hq = 2.7 and set a termination con- 
dition for the iteration, such as \nk+i - nk\ < e with e a small 
quantity. Once Eq. (6) is solved, the Hubble expansion rate 
can be obtained by 



H(z) = Ho 



^m0(l + Z? + ^r0(l + Z)^ + ^ko(l + zf 



1 - ^deiz) 



1/2 



(7) 

With the background evolution of the NADE model com- 
pleted, we now consider how to handle the perturbation of 
NADE. Since the EOS of NADE is always greater than -1, 
we can easily handle its perturbation. We take NADE as a per- 
fect fluid with the EOS given by Eq. (5). In the synchronous 
gauge, the conservation of energy-momentum tensor Ty.^ = 
gives the perturbation equations of density contrast and veloc- 
ity divergence for the NADE in the Fourier space. 



5' = -(l + Wde)(0+^)-3aH(c^,-Wde)S, (8) 



dz 



1+z 



(1 - ^de) 



3 -h G(z) ^|^de 



(6) 



where G(z) = 



^mO ( 1 +z)+flK) ( 1 +z)^+^kO 



with Q^o = ^dmo + 



^bO- 



ff = -aH{\-3c])e+ — k^S-k^cr. (9) 

1 + 

Here / = d/drj, H and Wde are given by Eqs. (7) and 
(5). Other notation is consistent with the work of Ma and 
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Bertschinger [12]. For the gauge ready formalism about the 
perturbation theory, see Ref. [13]. For the NADE, we assume 
the shear perturbation cr = 0. In our calculations, the adiabatic 
initial conditions will be taken. 

Finally, we feel that it may be necessary to make some ad- 
ditional comments on the treatment of dark energy perturba- 
tions in the NADE model. Actually, it is fairly difficult to 
calculate the cosmological perturbations of dark energy in the 
holographic-type dark energy models, because there is some 
non-local effect that makes the perturbation mechanism ex- 
tremely obscure and hard to handle. This is the reason why for 
a long time no one addressed the issue of dark energy pertur- 
bations in a complete manner in the holographic-type models. 
Be that as it may, Li, Lin, and Wang [14] made a calculation 
for the perturbations in the holographic dark energy model, 
in which some approximations are used in dealing with the 
non-local integrals in the equations governing the evolution 
of metric perturbations, and they proved that the perturbations 
are stable for both super-horizon and sub-horizon scales. But 
a complete analysis of the cosmological perturbations in the 
holographic dark energy model is still absent. Thus, in the 
past, to avoid the inclusion of dark energy perturbation prob- 
lem, one only used the CMB "distance priors" data, instead 
of the full CMB data, to constrain the holographic and age- 
graphic dark energy models. 

However, recently, it was realized that, in order to make 
progress, one should first ignore the non-local eff'ect in the 
holographic-type models and directly calculate dark energy 
perturbations in these models as if they are usual perfect flu- 
ids. Alone this line, global fit analyses on the holographic dark 
energy model were performed [8, 15]. It should be stressed, 
however, that such a treatment is only an expedient measure, 
in the case that we are not capable of completely handling 
the non-local eff'ect in the calculation of dark energy perturba- 
tions in these models. But, on the other hand, since the non- 
local eff'ect is not a dominant factor, it is believed that the fluid 
approximation in the holographic-type dark energy models is 
reasonable. 



III. GLOBAL FIT RESULTS 

In our global fit analysis, the parameter n of the NADE 
model is treated as a derived parameter, whose value can be 
obtained by the parameters Q^o, ^dmid and Ily^o via the Newton 
iteration algorithm, as mentioned in Sec. II. The "CosmoMC" 
package uses the physical densities of baryon ojiy = fltoh^ 
and cold dark matter codm = ^dmoh^ instead of flho and ^dmo- 
Thus, our most general parameter space vector is: 



TABLE I: Global fit results of the new agegraphic dark energy model. 



Class 



Parameter 



Best fit with errors 



Primary 



lOOQ^o/?^ 

T 



ris 

log[10iOA, 



uss 



032+0 



mu-m-m 



Derived n 

Ho (km/s/Mpc) 

^deO 
^mO 

Age (Gyr) 



n 1 ia+O.006+0.01 1+0.015 
^' ^ ^-^-0.005-0.009-0.014 

2 262'^^ -^^^'^^-'^^^'^^ ' 
0.020 
0.085; 

n Q7n+0-013+0;025+a036 
^•^^^-0.011-0.023-0.037 
o 1 oi +0.040+0.080+0.122 

z,.u / •^_o.077-0.15 1-0.222 
f.Q 9 + 1.2+2.6+3.9 
'J^-^-1.4-2.9-4.2 

697+0-011+0.025+0.037 



)14+0, 
„.J10-0.„__ 
+0.202+0.457+0.6 



0.284^^ 
13.125 

1 r\ C0 1 +i;238+2:485+3:696 
^^•^•^ ^-0.884-2.058-3.286 



The sound speed of dark energy, = Spde/^Pde, is usually 
treated as a phenomenological parameter for a fluid model of 
dark energy. It must be real and non-negative to avoid unphys- 
ical instabilities. In fact, if we treat the dark energy strictly as 
an adiabatic fluid, then the sound speed Cs would be imag- 
inary (in this case, the physical sound speed is equal to the 
adiabatic sound speed, = < 0), leading to instabilities 
in dark energy. In order to fix this problem, it is necessary to 
assume that dark energy is a non-adiabatic fluid and impose 

> by hand [16]. Our analysis is insensitive to the value 
of Cs. As long as Cs is close to 1, the dark energy does not 
cluster significantly on sub-Hubble scales. Therefore, in our 
analysis, we set Cs to be 1 ; the adiabatic sound speed can be 
imaginary, = dpde/dpde < 0. This is what is done in the 
CAMB and CMBFAST codes. Of course, one can also take 
Cs as a parameter, but the fit results would not be aff'ected by 
this treatment [17]. 

For the data, besides the CMB data (including the WMAP 
7-yr temperature and polarization power spectra [10]), the 
main other astrophysical results that we shall use in this paper 
for the joint cosmological analysis are the following distance- 
scale indicators: 

• The Union2.1 sample of 580 SNIa with systematic er- 
rors considered [18]. 

• The BAO data, including the measurements of 
rs(Zd)/Dv(0.2) and r,(zj)/Z)v(0.35) from the SDSS 
DR7 [19] and the measurements of A parameter at 
z = 0.44, 0.6, 0.73 from the WiggleZ Dark Energy Sur- 
vey [20]. 



V = (oJb, ojdm, O, T, Qko, ns. As), 



(10) 



where O is the ratio (multiplied by 100) of the sound horizon 
to the angular diameter distance at decoupling, r is the op- 
tical depth to re-ionization. As and ris are the amplitude and 
the spectral index of the primordial scalar perturbation power 
spectrum. For the pivot scale, we set kso = 0.002Mpc"^ to 
be consistent with the WMAP team [10]. Note that we have 
assumed purely adiabatic initial conditions. 



• The Hubble constant measurement Hq = 73.8 ± 
2.4km/s/Mpc from the WFC3 on the HST [21]. 

Our global fit results are summarized in Table I and Fig. 1. 
We notice two remarkable results: (i) the value of flko is 
greater than at least at the 3cr confidence levels (CLs), which 
implies that the NADE model distinctly favors an open uni- 
verse; (ii) the best-fit value of n is evidently less than those 
obtained in the previous work. 
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FIG. 1: The ID marginalized distributions of individual parameters and 2D contours at the l-3cr confidence levels from the global fit using 
the CMB+BAO+SNIa+//o data. 
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FIG. 2: Parameter spaces in the Qdeo-^mo plane constrained by dif- 
ferent joint data sets. The cyan, orange and red regions denote 
the constraint results from the joint CMB, BAO and Hq observa- 
tions without SNIa data, with Union2 sample of SNIa data and with 
Union2.1 sample of SNIa data, respectively. Note that the Union2.1 
sample is used with the systematic errors while the Union2 sample 
isn't, so it is no wonder why the Union2 sample slightly constrains 
the parameters more tightly than the Union2. 1 sample. 



The result of the spatial curvature we obtained is still at 
the level of 10"^, consistent with the other results about dark 
energy in the literature. However, it should be pointed out that 
our work is the first to show that the NADE model distinctly 
favors an open universe under current data. 



FIG. 3: The CMB Cj^ power spectra for the NADE and ACDM 
models with the best-fit parameter values. The black dots with error 
bars denote the observational data from the WMAP 7-yr observa- 
tions. 



9 ^7^+0.053+0.127+0.199 i_ 
z..u/j» 0.077-0.151-0.222 ^^^^ ^ 



The fit result of n is n 
3cr CLs. We find that the best-fit value of n is significantly 
diflTerent from the results obtained in the previous work [7, 
9, 22-25]. Eor example, Refs. [7], [22], [24] and [25] gave 
n = 2.810, 2.807, 2.887, and 2.886, respectively. What can 
this change tell us? Next, we shall explore what the factor is 
that is making this change. 

To do so, and without loss of generality, we take the most 
recent work in Ref. [7] as an example to make a comparison. 
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We list here four different points between them: (1) the spatial 
curvature is considered in this work while it isn't considered 
in Ref . [7] ; (2) the Hubble constant measurement and the 
BAO measurement of A parameter from the WiggleZ Dark 
Energy Survey are used in this work while they are not used 
in Ref. [7] ; (3) the SNIa data set used in this work is Union2. 1 
sample while it is Union2 sample in Ref. [7]; (4) this paper 
performs a global fit where the full CMB information from 
WMAP-7yr observations is used while the CMB information 
used in Ref. [7] comes only from 7-yr WMAP "distance pri- 
ors". After testing these four items one by one by varying 
one condition and setting the others the same, we find that the 
change is mainly due to the diff'erent SNIa data used. 

To see clearly, we plot the test results of diff'erent SNIa data 
used for the NADE model in Fig. 2. The cyan, orange and 
red regions denote the parameter spaces constrained by using 
the joint CMB, BAO and Hq observations without SNIa data, 
with Union2 sample of SNIa data, and with Union2.1 sample 
of SNIa data, respectively. Note that the CMB data used here 
are the 7-yr WMAP "distance priors" for simplicity. Com- 
paring the cyan region (without SNIa) with the orange region 
(with Union2), one can easily find that their overlapping re- 
gion is very small, which implies that there is a tension be- 
tween the Union2 sample of SNIa and the joint CMB, BAO 
and Hq observations. However, the tension between SNIa 
and CMB-hBAO-h//o will be greatly relieved once the Union2 
sample is replaced by the Union2.1 sample of SNIa, since the 
overlap between cyan and red regions is very large. Thus, 
we can conclude that the new constraint result of the NADE 
model obtained in this work is more reasonable than before. 

In addition, by using the global fit results obtained in this 
work (best-fit values), we plot the CMB power spectrum for 
the NADE model in Fig. 3. To make a comparison, we also 
plot the corresponding curve of the ACDM model with the 
best-fit parameters given by the same data sets. The 7-yr 
WMAP observational data points with uncertainties are also 
marked in the figure. We can clearly see that the NADE model 
produces a diff'erent feature compared with the ACDM model 
in the low / regions. However, since the low / data are rare and 
rough, the two models cannot be discriminated by the late- 
time ISW eff'ects by far. 



the diff'erential equation of Qj^, which can reduce one free pa- 
rameter for the NADE model. In order to easily modify the 
"CosmoMC" package, we set D^dmO, ^bo, and flko as free pa- 
rameters and ^ as a derived parameter, and use the Newton 
iteration algorithm to obtain the true value of n. 

Our global fit results show that the NADE model distinctly 
favors an open universe, since the value of flko is greater than 
at the 3(T CLs. Besides, we find that the value of the key 
parameter of the NADE model n = 2.673+^-^^^+^-}^^+0-i99 sig- 
nificantly diff'erent from those obtained in the previous work. 
We find that the main reason leading to such a change comes 
from the diff'erent SNIa samples used. Namely, the Union2.1 
sample and the Union2 sample give fairly diff'erent constraint 
results for the NADE model. Our further test indicates that 
there is a distinct tension between the Union2 sample of SNIa 
and other observations, and the tension will be removed once 
the Union2 sample is replaced by the Union2.1 sample. Thus, 
we conclude that the new constraint result of the NADE model 
obtained in this paper is more reasonable than the previous 
ones. 

We also plot the CMB power spectra for the NADE and 
ACDM models by using the best-fit results of the global fits. It 
is shown that the NADE model and the ACDM model produce 
slightly diff'erent features in the low / regions. Since the low 
/ data are rather rare and rough, the two models cannot be 
discriminated by the late-time ISW eff'ects yet. We expect that 
the future highly accurate data can do this job. 

If dark energy is truly dynamical, as hinted by recent 
work [26], the possible non-gravitational interaction between 
dark energy and dark matter may deserve deeper investiga- 
tions. In this work the interaction between dark energy and 
dark matter is not considered. If such an interaction exists in 
the NADE model, the calculations of the perturbations of dark 
energy and dark matter will become more complicated. The 
global fit analysis on the interacting agegraphic dark energy 
model is fairly attractive. In addition, the comparison of the 
various holographic models of dark energy under a uniform 
global fit is also necessary. We leave these subjects in the fu- 
ture work. 



IV. CONCLUDING REMARKS 

In this work, our main task is to perform a global fit anal- 
ysis on the NADE model in a non-flat universe by using a 
MCMC method with the joint data of the full CMB spectra, 
BAO, SNIa and Hq observations. We use the initial condition, 

^deizini) = ^^^^i^ (l + V^Xi^)' at ZM = 2000, to solve 
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